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Abstract 

The treatment of anharmonic oscillators (including double wells) by instan- 
ton methods is wellknown. The alternative differential equation method - 
excepting various WKB eigenvalue approximations — is not so wellknown. 
Here we reformulate the latter completely parallel to the strong coupling case 
of the cosine potential and Mathieu equation for which extensive literature 
and monographs exist. The solutions and eigenvalues of the anharmonic os- 
cillator equations are obtained analogously as asymptotic expansions, and 
the exponentially small deviations are shown to follow from the imposition 
of boundary conditions — completely parallel to the case of the Mathieu 
equation. In spite of the additional involvement of WKB solutions, the ex- 
ploitation of parameter symmetries of the basic equation (as in the case of 
the Mathieu equation) leads to a double well level splitting which agrees with 
the otherwise Furry-factor improved WKB approximation (and, of course, 
with calculations using periodic instantons). 

1 Introduction 

In a recent paper on new ways to solve the Schrodinger equation, R. Friedberg 
and T. D. Lee P described the anharmonic oscillator "long stand- 
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ing difficult problem of a quartic potential with symmetric minima" . Quartic 
anharmonic oscillator potentials are for many reasons of fundamental inter- 
est in quantum mechanics and have therefore repeatedly been the subject 
of detailed investigations in diverse directions. In particular the investiga- 
tions of Bender and Wu [2], [H| which related analyticity considerations to 
perturbation theory and hence to the large-order behaviour of the eigenvalue 
expansion attracted widespread interest. An important part of their work 
was concerned with the calculation of the imaginary part of the eigenen- 
ergy in the non-self adjoint case which permits tunneling. Derivations of such 
a quantity are much less familiar than calculations of discrete bound state 
eigenenergies in quantum mechanical problems. This lack of popularity of 
the calculation of complex eigenvalues even in texts on quantum mechanics 
may be attributed to the necessity of matching various branches of eigen- 
functions in domains of overlap and to the necessary imposition of suitable 
boundary conditions, both of which make the calculation more difficult. 

The development and recognition of the significance of the path-integral 
method, particularly in connection with consideration of pseudo-particle con- 
figurations, displaced the Schrodinger equation into the shadow of the path- 
integral and suggested that little else would be gained from further studies 
of the differential equation. In addition, the equivalence of the two methods 
(at least in the one loop approximation) is not always appreciated [I]. 

Our intention in the following is to present treatment of the three different 
types of quartic anharmonic oscillators (these types arising from different 
signs of quadratic and quartic terms) formulated in parallel to the strong 
coupling case of the cosine potential 0,0, the analogy with the periodicity 
of the latter entering only in the formulation of the boundary conditions. 
Indeed the Schrodinger equation with quartic oscillator potential should be 
raised to the rank of a standard differential equation akin to the Mathieu 
and similar equations which are not reducible to some hypergeometric type. 
Thus compare the equations 



Large-/z 2 asymptotic expansions of the Mathieu equation are wellknown, and 
their lifting of the degeneracy of the asymptotically degenerate oscillator 
levels (i.e. for h 2 — > oo) yields the boundaries of the energy bands obtained 




y"(z) + [E- V(z)]y(z) = 
J 2h 2 cos 2z = 2h 2 - Ah 2 z 2 + §/iV • ■ ■ 



Mathieu, 
quartic. 
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with the (convergent) small- h 2 expansion [Zj,[Hl- Books on Special Functions 
of Mathematical Physics contain all sorts of asymptotic expansions like the 
large-/? 2 expansion above, so that there is no reason to see in the allegedly 
"divergent perturbation series" of the anharmonic oscillator an ill-natured 
problem. These series are completely analogous to asymptotic expansions of 
Special Functions. 

The treatment of the double- well potential by instanton methods and the 
calculation of its lowest level splitting jU] constitute a standard example of 
the path-integral method using instantons. The corresponding calculation 
for higher levels requires a different method ^0], [H] |JZ| and is achieved 
with the use of periodic instantons o r other multi-instanton methods 
|14j . For the inverted double- well potential corresponding path- integral cal- 
culations using periodic instantons have been carried out in Ref. j!5j . It 
is natural therefore to compare the Schrodinger and path-integral methods 
and it is reassuring to obtain in both cases identical results. Earlier WKB 
investigations of the eigenvalues of the Schrodinger equation with double well 
potential serve as valuable checks of the results of other methods, particularly 
since these have in some cases also been compared with numerical methods. 
However pure WKB results do not always agree with one-loop path-integral 
results jn]. Thus improved WKB approximations to double well eigenvalues 
have in particular been investigated and results in agreement with the instan- 
ton method have been achieved. Such WKB method investigations can be 
found particularly in Refs. [TS], [Hj, [IE], [JH] (regrettably with considerable 
suppression of intermediate steps). 

In the following we review the method of the differential equation |2Uj . 
however reformulated in parallel to the corresponding treatment of the Math- 
ieu equation [3] for which the perturbation method in terms of a tunneling 
deviation from an integer (in the following q — qo, go = 1, 3, 5, . . . ,) was origi- 
nally developed. In our treatment here we leave various detailed calculational 
aspects to appendices. A feature different from that of the Mathieu equa- 
tion is the additional involvement of WKB solutions for the evaluation of 
the boundary conditions at the central maximum in the case of the double 
well potential. However, exploiting parametric symmetries of the original 
Schrodinger equation, the level splitting obtained is in complete agreement 
with the Furry-factor corrected WKB result and that of the one-loop approx- 
imation in path-integral methods, the Furry-factor being a correction factor 
to the normalisation constant of the WKB wave function [21]. More detailed 
considerations are collected in our appendices which contain in addition to 



3 



calculational aspects, some WKB formulae for easy reference, since these are 
not always available but are are essential in this context and are frequently 
or casually referred to as easily obtainable. 



2 The Three Types of Potenials 

In the cases treated most frequently in the literature the anharmonic oscilla- 
tor potential is defined by the sum of an harmonic oscillator potential and a 
quartic contribution. These contributions may be given different signs, and 
thus lead to very different physical situations, which are nonetheless linked 
as a consequence of their common origin which is for all one and the same 
basic differential equation. 




(1) (2) (3) 

Fig. 1 The three different types of anharmonic potentials. 

To avoid confusion we specify first the potential V(z) in the Schrodinger 
equation 

^ + [E~ V(z)]y(z) = 

for the different cases which are possible and illustrate these in Fig. 1. The 
three different cases are: 

(1) Discrete eigenvalues with no tunneling: In this case 

V ( z ) = -\h A \z 2 + -\c 2 \z 4 , 

(2) Discrete eigenvalues with tunneling: In this case, described as the case 
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of the double well potential, 



V(z) = -hh i \^+-\d i \z 4 - 

w 4 2 



(3) Complex eigenvalues with tunneling: In this case, with the potential 
described as an inverted double well potential, 

V {z) = -\h 4 \z 2 --\c 2 \z\ 

Case (1) is obviously the simplest with the anharmonic term implying 
simply a shift of the discrete harmonic oscillator eigenvalues with similarly 
normalisable wave functions. The shift of the eigenvalues is easily calcu- 
lated with straightforward perturbation theory. The result is an expansion 
in descending powers of h 2 . It is this expansion which led to a large number 
of investigations culminating (so to speak) in the work of Bender and Wu 
|2],jSI w ho established the asymptotic nature of the expansion* and its Borel 
summability.t The resulting eigenvalues are given by Eq. (34) below with 
Q — Qo — 2n + 1, n = 0, 1, 2, . . . and c 2 replaced by — c 2 , i.e. 

E(q, h 2 ) = \ % h 2 + ^(g 2 + 1) - |^(4g 3 + 29g ) + O . 

Case (2) is also seen to allow only discrete eigenvalues (the potential 
rising to infinity on either side), however the central hump with troughs on 
either side permits tunneling and hence (if the hump is sufficiently high) a 
splitting of the asymptotically degenerate eigenvalues in the wells on either 
side which vanishes in the limit of an infinitely high central hump. The 
resulting eigenvalues are given by Eq. (145) below, i.e. 

2g +i/ i 2( JilW2 h e 
E( go , h 2 ) . E ( go , h 2 ) =F ^^[Ifa ,!)]! 6 "^' * = 1, 3, 5, ... , 

where E (qo, h 2 ) is given by Eq. (}95|) or Eq. (|C.1J) . i.e. 

F(a h 2 )- ^ + 1 ah 2 C2(3g ° 2+1) I Of 1 



*This is shown by demonstrating that the i-th late term behaves like a factorial in % divided by an 
i-th power. See Ref. 1221 . 

^This is the case when the late terms alternate in sign. 
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Case (3) is seen to be very different from the first two cases, since the 
potential decreases without limit on either side of the centre. The boundary 
conditions are non-selfadjoint and hence the eigenvalues are complex. This 
type of potential allows tunneling through the barriers and hence a passage 
out to infinity so that a current can be defined. If the barriers are sufficiently 
high we expect the states in the trough to approximate those of an harmonic 
oscillator, however with decay as a consequence of tunneling. The resulting 
complex eigenvalues are given by Eq. (76) together with Eq. (34), i.e. 



The question is therefore: How does one calculate the eigenvalues in these 
cases with the help of the Schrodinger equation? This is the questions we 
address here, and we present a fairly complete treatment of the case with large 
values of h 2 along lines parallel to those in the case of the cosine potential in 
Ref. We do not dwell on Case (1) since this is effectively included in the 
first part of Case (3), except for a change of sign of |c 2 |. Thus we are mainly 
concerned with the double well potential and its inverted form.* We begin 
with the second. In this case our aim is to obtain the aforementioned complex 
eigenvalue. In the case of the double well potential our aim is to obtain 
the separation of harmonic oscillator eigenvalues as a result of tunneling 
between the two wells. Calculations of complex eigenvalues (imaginary parts 
of eigenenergies) are rare in texts on quantum mechanics. We therefore 
consider here in detail a prominent example and in such a way, that the 
general applicability of the method becomes evident. 

'''In each of these two cases we have to take two different boundary conditions into account. In the 
cosine case of Ref. only one type is required at a = w/2. The reason is that this is actually derived 
from the boundary condition at z = using the periodicity of this case; see Ref. 0, equations (6), p. 
108, and the relations (9), p. 100. 




with 



Eo(qo, h 2 ) = l -q Q h 2 - ^(q 2 + 1) - -^(4g 3 + 29q ) + O 
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3 The Inverted Double- Well Potential 



3.1 Defining the problem 

We consider the case of the inverted double- well potential depicted as Case 
(3) in Fig. 1. The potential in this case is given by 

V(z) = -v(z), v{z) = -\h i z 2 + l -c 2 z\ (1) 

for h 4 and c 2 real and positive, and the Schrodinger equation to be considered 
is 

^ + [E + v(z)]y = 0. (2) 




Fig. 2 The inverted double-well potential with (hatched) oscillator potential. 

We adopt the following conventions which it is essential to state in order 
to assist comparison with other literature. We take h = 1 and the mass 
mo of the particle = 1/2. This implies that results for mo = 1 (a frequent 
convention in field theory considerations) differ from those obtained here by 
factors of 2 1 / 2 , a point which has to be kept in mind in comparisons. If 
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suffixes 1/2, 1 refer to the two cases, we can pass from one to the other by 
making the replacements: 

E 1/2 = 2Ei, hf /2 = 2hi, c 2 /2 = 2c 2 . (3) 

Introducing a parameter q and a quantity A = A(q,h), and a variable w 
defined by setting 

E = ^qh 2 + — and w = hz, (4) 
we can rewrite Eq. (J2J) as 

V g (w)y(w) = -1(A + c 2 w 4 )y(w) (5) 

with 

In the domain of w finite, \h 2 \ — > oo and c 2 finite, the harmonic part of the 
potential dominates over the quartic contribution and Eq. (jSJ) becomes 

V q (w)y(w) = o(-^j. (7a) 

The problem then reduces to that of the pure harmonic oscillator with y(w) 
a parabolic cylinder function, i.e. 

y(w) oc Dx^^w) and q — q — 2n+l, n — 0, 1, 2, (7b) 

The perturbation expansion in descending powers of fa suggested by the above 
considerations is therefore an expansion around the central minimum of V(z) 
at z — 0. The positions z± of the maxima of V(z) on either side of z = in 
the case c 2 > are obtained from 

h 2 

v'(z±) = as z± — ±— (8) 

with 

fa 8 

i/'(*±) = fa 4 and V(*t) = 

2°er 
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Thus for c 2 > and relatively small, and h 2 large the eigenvalues are es- 
sentially perturbatively shifted eigenvalues of the harmonic oscillator as is 
evident from Fig. 1. 

The problem here is to obtain the solutions in various domains of the 
variable, to match these in domains of overlap, then to specify the necessary 
boundary conditions and finally to exploit the latter for the derivation of the 
complex eigenvalue. The result will be that derived originally by Bender and 
Wu 05 jS], although our method here (which parallels that used in the case 
of the cosine potential) is different. 

3.2 Three pairs of solutions 

We are concerned with the equation 



d 2 y{z) 
dz 2 



+ 



h 4 z 2 
E- — + 



where 



E 



2 

A 



-qh 2 + - 
T 2h 4 



(9) 



(10) 



Here again q is a parameter still to be determined from boundary conditions, 
and A = A(g, h) is obtained from the perturbation expansion of the eigen- 
value, as encountered and explained earlier. Inserting (jlOj) into we obtain 



dz 2 



+ 



1 , 2 A h^z 2 
2 qk + 2i?- — + 



c 2 z 4 



y = 0. 



(11) 



The solutions in terms of parabolic cylinder functions are valid around z = 
and extend up to z ~ 0(l/h 2 ), as we shall see. Before we return to these 
solutions we derive a new pair which is valid in the adjoining domains. Thus 
these solutions are not valid around z = 0. In order to arrive at these 
solutions we set in Eq. (fTTj) 

l/2-i 



y{z) = A(z) exp 



±i 



dz 



h A z 2 



c 2 z A 



4 2 

Then A(z) is found to satisfy the following equation 



A"(z) ± 2i 



h A z 2 



c 2 z 4 



1/2 



A\z)±iA(z)^- 



4^2 



h 4 Z 



c 2 z* 



+ 



1 L 2 A 

-qh + — tt 

2 H 2h A 



A(z) 



0. 



'121 



1/2 



(13) 
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Later we will be interested in the construction of wave functions which are 
even or odd around z — 0. This construction is simplified by the consider- 
ation of symmetry properties of our solutions which arise at this point. We 
observe — before touching the square roots in Eq. (|T3*|) — that one equation 
(of the two alternatives) follows from the other by changing the sign of z 
throughout. This observation allows us to define the pair of solutions 



y A (z) = A(z)exp 
Va( z ) = M z ) ex P 



h A z 
- + 

h 4 z 2 



C 2 Z 4 1/2-, 



2 4 ^1 1/2- 
C Z Z* ' 



(14a) 

dz{ - — + j, (14b) 

A(z) = A(-z) and y A (z) = y A {-z), (15) 

where A(z) is the solution of the upper of Eqs. (fTSj) and A(z) that of the 
lower of these equations. We take the square root by setting 



with 



zW ^ 1/2 



ztf 



For large h we can write the Eqs. (|13|) as 

Tz A'(z) T \A(z) + ^qA(z) = o(^- 2 
We define A q (z) as the solution of the equation 

zA' q (z)-^(q-l)A q (z) = 0, 



(16) 



(17) 



i.e. 



A q (z) 



zz 



1(9-1) 



( 2 2)l/4 



exp 



dz 



^2)1/2 



We define correspondingly 

A g (z) = z'^+V = A_ q {z) 



exp 



-Q 



dz 



^2)1/2 



(19) 



We see that one solution follows from the other by replacing z by —z. Clearly 
A q (z),A q (z) approximate the solutions of Eqs. (fT3|) and we can develop a 
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perturbation theory along the lines of our method as employed in the case 
of periodic potentials. One finds that these solutions are associated with 
the same asymptotic expansion for A — to be derived in detail below in 
Appendix A and as a verification again in connection with the solution ub 
- as the other solutions. Since these higher order contributions are of little 
interest for our present considerations, we do not pursue their calculation. 
Thus we now have the pair of solutions 



y A (z) = exp 



dz 



1 2 1,4 , 1 2 4 
-z h + -C Z 

4 2 



V2i r 



A q {z) + 



h 2 



(20a) 



VaK z ) 



exp 



dz 



1 2t4 ,-'-2 4 

—z h H — c z 
4 2 



1/2' 



A q (z) 



O 



. (20b) 



These expansions are valid as decreasing asymptotic expansions in the do- 
main 

1 



i.e. away from the central minimum. With proper care in selecting signs of 
square roots we can use the solutions ()20a|) and ()20b|) to construct solutions 
y±(z) which are respectively even and odd under the parity transformation 



z — > —z (or equivalently q — > —q, h 2 

1 



-h ), i.e. we write 



y±\ z ) 



[y A (q,h ;z) ±y A (q,h ; z)}. 



(21) 



In the next two pairs of solutions the exponential factor of the above 
solutions of type A is contained in the parabolic cylinder functions (which 
are effectively exponentials times Hermite functions). 

We return to Eq. (jSJ). The solutions y q (z) of the equation 



V q (w)y q (w) =0, w = hz, 



(22) 



are parabolic cylinder functions Di^ q _^(±w) and D_i^ q+1 -^(±iw) (observe 
that Eq. ()22|) is invariant under the combined substitutions q — > —q, w — * 
±iw) or functions 



B q (w) 



[l( g -l)]!2ifo-D 



and C q (w) 



D_ l{q+1) {±tw)2^ 



(23) 
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These solutions satisfy the following recurrence relation (obtained from the 
basic recurrence relation for parabolic cylinder functions given in the litera- 
ture^) 

w 2 y q = -(q + 3)^+4 + qy q + -(q~ 3)y ff _ 4 . (24) 

The extra factors in Eq. (|23|) have been inserted to make this recurrence 
relation assume this particularly symmetric and appealing form.^ For higher 
even powers of w we write 



w 2i 



y q =Yl s ^(<iAj)y q +4j, (25) 



j=-l 

where in the case % — 2 

S 4 (q,±8) = i(g±3)(g±7), 
S,(q,±4) = (q±2){q±3), 

Si(g,o) = ^ 2 + i)- (26) 

The first approximation y(w) = y(°\w) = y q (w) therefore leaves uncompen- 
sated terms amounting to 

= -i(A + c 2 w*)y q (w) = -1 £ [g,g + 4j]y ?+4i , (27) 

i=-2 

where 

[q,q] = A + (?S i (q,0), andforj^O: [q, q + 4j] = c 2 S 4 (g, 4j). (28) 



^Ref. |23l . pp. 115 - 123. Comparison with our notation is easier if this reference is used, 
an alternative to B q (w) in Eq. 1231 one can choose the solutions as B q (w) with 

B q (w) 



[i(g-3)]!24(9-D 
These satisfy the recurrence relation 



u 2 y q (w) = ~(q + l)y q+ A + qy q + ~(<J - l)y q -4- 
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Now, since V q y q = 0, we also have V q+ 4jy q+4: j = 0, and so 

V q y q+ Aj{w) = -Ajy q+4j (w). (29) 

Hence a term ^y q +Aj on the right hand side of Eq. (|2*7jl can be removed by 
adding to the contribution (— ^/4j)y g+4 j. Thus the next order contribu- 
tion to y q is 

y () ( w ) = w 1^ — t- — Vi+V- ( 3 °) 

For the sum = y^°\w) + y^\w) to be a solution to that order we must 
also have to that order 

[q, q] = 0, i.e. A = ^-{q 2 + l)c 2 + O f i) . (31) 
Proceeding in this way we obtain the solution 

y = yW) ( W ) + yW ( W ) + y( 2 ) («;) + ... (32) 

with the corresponding equation from which A can be obtained, i.e. 

= i[ 9 , 9 ] + (^) 2 g^^[ 9 + 4i, 9] + .... (33) 

Evaluating this expansion and inserting the result for A into Eq. (fTUJl we 
obtain 

E i«> k ") = \^ ~W iq2 + 1) ~ + ^ + (^6 ) • (34) 

We observe that odd powers of q arise in combination with odd powers of 
1/h 2 , and even powers of q in combination with even powers of 1/h 2 , so that 
the entire expansion is invariant under the interchanges 

q — > —q, h — »■ — h 2 . 

This type of invariance is a property of a very large class of eigenvalue prob- 
lems. Equation (|34jl is the expansion of the eigenenergies E of Case (1) 
with q — q — 2n + 1, n — 0, 1, 2, . . . . In Case (3) the parameter q is only 
approximately an odd integer, the difference q — qo arising from tunneling.. 
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We can now write the solution y(w) in the form 



y( w ) = y q i w ) + 



oo / j \ i 2i 
i=l ^ / j=-2i,j^0 



W 



(35) 



where for instance 
Pi{q,q±4) 

P2(q,q±4:) 



[g,g±4] ^ c2 (g±2)(g±3) 



±4 



±4 



[g,g±4] [g±4,g±4] + [g,g±8] [g±8,g±4] 



±4 ±4 
[g, g =f 4] [g =f 4, g ± 4] 



±4 



T4 



±4 



and so on. Again we can write down a recurrence relation for the coefficients 
Pi(q,q + 4j), i.e. 



4tPi(q, q + 4t)=J2 9 + 4j + 4t)[g + 4j + 4t, g + 4t] (36) 

j=-2 

with the boundary conditions 

Po(q,q) = 1) and for j ^ all other Po(q, q + 4j) = 0, 
PiM<li l) = °' 

Pi(q, g + 4j) = for |j| > 2z or |j| > 2z + 1. (37) 

For further details concerning these coefficients, their recurrence relations and 
the solutions of the latter we refer to Ref. [20] • Since our starting equation 
(II lj) is invariant under a change of sign of z, we may infer that given one 
solution y(z), there is another solution y(—z). We thus have the following 
pair of solutions 



Vb{z) 



yB\ z ) 



co j ^ v i 2i 

+ E Pi(q,q + 4j)B g+ *jH 

i=l ^ ' ,=-2i,7^0 



[|/B(2)]arg^=7r = fe/fl(-z)]arg*=0- 



ui=/iz,arg 2=0 

(38 



These solutions are suitable in the sense of decreasing asymptotic expansions 
in the domains 
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They are linearly independent there as long as q is not an integer. 

Our third pair of solutions is obtained from the parabolic cylinder func- 
tions of complex argument. We observed earlier that these are obtained by 
making the replacements 

q — > — q, w — » ±iw. 
These solutions are therefore defined by the following substitutions: 

yc(z) = [y B (z)] q ^- q , h ^ih, y c ( z ) = [y B ( z )}q^~q,h^ih (39) 

with the same coefficients Pi(q,q + 4j) as in y B . The solutions yciVc are 
suitable asymptotically decreasing expansions in one of the domains 

1 \ 71 

We emphasise again that all three pairs of solutions are associated with 
the same expansion of the eigenvalue E(q, h 2 ) in which odd powers of q are 
associated with odd powers of h 2 , so that the eigenvalue expansion remains 
unaffected by the interchanges q — > —q, h 2 — > —h 2 as long as corrections 
resulting from boundary conditions are ignored. 

3.3 Matching of solutions 

We saw that the solutions of types B and C are valid around the central 
minimum at \z\ = 0, the solutions of type A being valid away from the 
minimum. Thus in the transition region some become proportional. In order 
to be able to extract the proportionality factor between two solutions, one has 
to stretch each by appropriate expansion to the limit of its domain of validity. 
In this bordering domain the adjoining branches of the overall solution then 
differ by a constant." 

First we deal with the exponential factor in the solutions of type A. Since 
these are not valid around z = but h 2 is assumed to be large, we expand 

^Variables like those we use here for expansion about the minimum of a potential (e.g. like w of 
Eq. 0) are known in some mathematical literature as "stretching variables" and are there discussed in 
connection with matching principles, see e.g. Ref. 1241 . 



<o 
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the integrand as follows: 
exp 

= exp 

= exp 
= exp 



1 



dz{ - =-z 2 h A + '-c A z 



1 



,2 A 



4 2 
2cV 



1/2- 



h^2 
8c 2 3 

12c 2 



h 4 



2c 2 z 2 ^ 3/2 



1 - 



2cV ^ 1/2 1 



Z; 1 



(40) 



Considering the pair of solutions Ua{z), Va( z ) we see that in the direction of 
z = (of course, not around that point) 



Va{z) 
Va{z) 



3 ft 6 /12c : 



e 4- 



e _ fc 6 /12c 2 e l z 2 h 2 



2-1(3+1)+ Of 1 



(41) 



The cases of the solutions of types 5 and C require a careful look at 
the parabolic cylinder functions since these differ in different regions of the 
argument of the variable z. Thus from the literature j2H] we obtain 



i=0 Z 'L2 



y [|(g -l)]» 



but 



,-iri 



4z-l)]!(-2w 2 ) i ' 



argw| < -7r, 



(42) 



i=0 



m(q-M-l)]\ (-2w 2 Y 



1(9 + 1)]! 



(2 7 r) 1 / 2 e- i f^- 1 ) 
[-1(9 + 1)]! W *(*+D ^ i!h|( ? + 4z + 1)]! (2w 



■ , 5 1 
with -7T > argw; > -7r. 

4 4 

The function Di ( ^(w) has a similarly complicated expansion for 



(43) 



1 5 

> argw; > — ir. 

4vr & 4 



1(3 



From 



we obtain for the solution ys(z),w = hz: 
(/i¥)iH) 



Vb{z) 



B q {w) 



[I( ? -1)]! 2 K^) 



-±h 2 z 2 

e 4 



(44) 



In the solution y B (z), with z in ?/b(z) replaced by — z, we would have to 
substitute correspondingly the expression ()43|) (since z — » —2; implies argz = 
±7r). We do not require this at present. Comparing the solution Ua{ z ) of 
Eq. (|4ip with the solution y B {z) of Eq. (|44j) we see that in their common 
domain of validity 

Va(z) = -y B (z) (45) 



with 



a 



2^(9-1) 



a 



e i2c^ 



1 + 



(46) 



[i(g-l)]!2^) 

However, the ratio of y A (z),y B (z) is not a constant. 

We proceed similarly with the solutions yc(z),y c (z). Inserting the ex- 
pansion ()42|) into yd z ) we obtain 



Vc\ z ) 



(_ /i 2 z 2)-I( g +l) e |^^ 
[_l(g+l)]!2-i(»+l) 



(47) 



Comparing this behaviour of the solution y c (z) with that of solution y A {z) 
of Eq. (|4~Tf , we see that in their common domain of validity 



a 



where 



a 



+ l)]!2-3^ +1 ) . 



1 + 



(48) 



(49) 



Again there is no such simple relation between ^(z) and yc(z) 



3.4 Boundary conditions at the origin 

(A) Formulation of the boundary conditions 

The really difficult part of the problem is to recognise the boundary con- 
ditions we have to impose. Looking at the potential we are considering here 
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- as depicted in Fig. 2 — we see that near the origin the potential behaves 
like that of the harmonic oscillator, in fact, our large-/? 2 solutions require this 
for large h 2 . Thus the boundary conditions to be imposed there are the same 
as in the case of the harmonic oscillator for alternately even and odd wave 
functions. Recalling the solutions y±(z) which we defined with Eq. (121)1 as 
even and odd about z = 0, we see that at the origin we have to demand the 
conditions 

y ' + (0) = and y_(0) = (50) 

and y+(0) ^ 0, y'_(0) ^ 0. The first of the conditions (|5D|) will be seen to 
imply q = 2n + 1 = 1, 5, 9, . . . and the second q = 3, 7, 11, . . .. For instance 
go = 1 (or = 0) implies a ground state wave function with the shape of a 
Gauss curve above z — 0, i.e. large probability for the particle to be found 
thereabouts. At z = the solutions of type A are invalid; hence we have to 
use the proportionalities just derived in order to match these to the solutions 
valid around the origin. Then imposing the above boundary conditions we 
obtain 



y'+(0) = lim-[y' A (z) +y A (z)] 

z— >0 2 



y B (o)+y c (o) 

a a 



(51) 



and 



y_(0) = \im-[y A (z) 

z— >0 2 

5fi»»(0)-^ o (0) 



(52) 



Thus we obtain the equations 



— = and 

a 



Vb(0) 
Vo(0) 



a 
a 



(53) 



Clearly we now have to evaluate the solutions involved and their derivatives 
at the origin. We leave details to Appendix B. 

(B) Evaluation of the boundary conditions 



We now evaluate Eqs. (|53j) in dominant order and insert from Eqs. 
P9~j) the appropriate expressions for a and a. Starting with the derivative 
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expression we obtain (apart from contributions of order 1/h 2 ) 

lsin{f(g + 3)} (/ i 2 )3("- 1 )[-i(g + l)]!2-3^+ 1 ) _*fi 



i sin{f (g - 3)} " [l( q - l)]!2a(«-i)(-/ l 2)-|(9+i) 



e 6 C ^ 



L 1 

We rewrite the left hand side as 

1 sin{f(g + 3)} 



7T 



i tan< — (g + 3) 



i sin{f (g + 3 -6)} [4 

We rewrite the right hand side of the derivative equation again with the help 
of the inversion and duplication formulae and obtain 



-7T- 



- 1)1 ! S(9 - 3)]! sin{f (g +1)} V 2 [|(g - 1)]! cos{f (q + 3)} 

Then the derivative relation of (|53|) becomes 

fvr, ,1 /F(/i 4 )^ 4 (-l)i(9+i) h6 

8in |_ (? + 3 )} = -^-LJ i _i_i_ e ^. ,54, 

Proceeding similarly with the second of relations (jf)3|) . we obtain 

cos |_ (9 + 3 )| = -^LJ_i_J_ e ^. (55) 

In each of Eqs. (J54)) and (|55|) the right hand side is an exponentially small 
quantity. In fact the left hand side of (J53|) vanishes for q — q — 1, 5, 9, . . . and 
the left hand side of (J53j) for q — go — 3, 7, 11, . . .. With a Taylor expansion 
about go the left hand side of ({MJ) becomes 

(? - 9b)J cos {f (?o + 3)} + • • • ~ (g - g )^(-ir^°+ 3 ). 
It follows that we obtain for the even function with g = g = 1, 5, 9, . . . 

2^2 (/l 2 W 2 . . 

g - g ~ ± J ni , e'g. 56 

A [5(90 - I)]' 

Expanding similarly the left hand side of Eq. f|55[) about go = 3, 7, 11, . . ., we 
again obtain ()56|) but now for the odd function with these values of go- We 
have thus obtained the conditions resulting from the boundary conditions at 
z = 0. Our next task is to extend the solution all the way to the region 
beyond the shoulders of the inverted double well potential and to impose 
the necessary boundary conditions there. Thus we have to determine these 
conditions first. 
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3.5 Boundary conditions at infinity 

(A) Formulation of the boundary conditions 

We explore first the conditions we have to impose at \z\ — > oo. Recall 
the original Schrodinger equation (j2J) with potential For c 2 < and the 
solution y(z) square integrable in — oo < !Rz < oo, the energy E is real; this 
is the case of the purely discrete spectrum (the differential operator being 
selfadjoint for the appropriate boundary conditions, i.e. the vanishing of 
the wave functions at infinity). This is Case (1) of Fig. 1. The analytic 
continuation of one case to the other is accomplished by replacing ±c 2 by 
=pc 2 or, equivalently, by the rotations 

E -> e iw E = -E, z-* e i7r/2 z, z 2 -> -z 2 . 

One can therefore retain c 2 as it is and perform these rotations. It is then 
necessary to insure that when one rotates to the case of the purely discrete 
spectrum without tunneling, the resulting wave functions vanish at infinity 
and thus are square integrable. Thus in our case here the behaviour of the 
solutions at infinity has to be chosen such that this condition is satisfied. 
Now, for — > ±oo we have 



y(z) ~ exp < ± % \ dz 



r 2 ' 

4 

— Z 

2 



1/2 ^ f / C 2X l/2 3 



In order to decide which solution or combination of solutions is compatible 
with the square integrability in the rotated (c 2 reversed) case, we set 



z = \z\e 



Then 



c 2 \ l/2 z n r /c 2 \ 1/2 h 3 



GXP ^ ~ H 2 J = 6XP \- 1 2 — 6 



3^6»^ 



fc 2 \ 1/2 \z\ 3 
exp { — il — J -^-(cos 36 + % sin 3$ 

c 2 \ 1/2 M 3 1 

K ex P \ ( ~2 j ~ sin 36 f ■ 
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c 2N 


f 


z 3 


2, 




y 


c 2 ^ 


f 




2 , 







This expression vanishes for \z\ — * +00 if the angle # lies in the range — tt < 
39 < 0, i.e. if -tt/3 < < 0. Thus 

exp ^ — « ( 77 ) — )■ ^ for JRz — > +00 in arg zG ^ — — , , , 

exp|+z^— J — j> ^0 for JRz — > —00 in argz G ^0,—^. (58) 

Rotating z by n/2, i.e. replacing sin 39 by 

sin 3^0 + I J = - cos 39 

we see that the solution with the exponential factor is exponentially decreas- 
ing for \z\ — > 00 provided that 

cos 39 > 0, 
i.e. in the domain — 7r/2 < 39 < tt/2, or 



TX TX 

— < 9 < -. 
6 6 



In the case of the inverted double well potential under consideration here 
(i.e. the case of complex E), we therefore demand that for Jiz — > +00 and 
—71/ 3 < argz < the wave functions have decreasing phase, i.e. 



c 2 \ 1/2 z 3 



^)~ex P j-*l-j -j, c 2 >0. (59) 

This is the boundary condition also used by Bender and Wu |2j- For c 2 < 
we have correspondingly 

y(z) ~ exp < ± I -— I — > , c 2 < 0, for z — > ^00. 



This is not the asymptotic behaviour of a wave function of the simple har- 
monic oscillator. We have to remember that we have various branches of the 
solutions y(z) in different domains of z. 

Our procedure now is to continue (in the sense of matched asymptotic 
expansions) the even and odd solutions (|2*Tj) to + infinity and to demand that 
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they satisfy the condition (|59|) for c 2 > 0. Equating to zero the coefficient of 
the term with sign opposite to that in the exponential of Eq. (|59|) will lead to 
our second condition which together with the first obtained from boundary 
conditions at the origin determines the imaginary part of the eigenvalue E. 



V(z) 





E >V 


/ E <V\ 




real 


/ real \ 










O z Z \ 
u 1 \ 



z 



Fig. 3 The inverted double-well potential with turning points zq,zi. 

(B) Evaluation of the boundary conditions 

The following considerations (usually for real z) require some algebraic 
steps which could obscure the basic procedure. We therefore explain our 
procedure first. Our even and odd solutions y±(z) (cf. Eq. (j2"T]) ) were defined 
in terms of solutions of type A which have a wide domain of validity. Looking 
at Fig. 3 we see that at a given energy E and to the right of z = (which is the 
only region we consider for reasons of symmetry) there are two turning points 
Zo, Z\. Thus we have to match the solutions of type A first to solutions to the 
left of Z\ and then extend these to solutions to the right, and there impose 
the boundary condition (JoTJj) on y±{z) (by demanding that the coefficient of 
the solutions with other behaviour be zero). We do this extension with the 
help of WKB solutions, i.e. we match the WKB solutions to the left of z\ 
to the solutions of type A, and then use the WKB procedure (called "linear 
matching" across the turning point) to obtain the dominant WKB solutions 
beyond Z\. 

The distant turning point at Z\ as indicated in Fig. 3 is given by (using 
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Eq. for E and ignoring nondominant terms) 



1 2t4 i 1 2 4 

•— 2 ft |-cz 
4 2 



i.e. 



Zi 



1 



2gc 5 



(60) 



The WKB solutions have been discussed in the literature.** We obtain in 
the domain V > dtE to the left of Z\ as the dominant terms of the WKB 
solutions 



Vwkb( z ) 



x exp 



x exp 



i 2 i 2t4 1 2 4 

-on. H — 2; n c z 

2 y 4 2 



-1/4 



dz 



l qh * + 1^4 _ 1 2 4 
2^ 4 2 



nl/2 



i o/i 2 + 1^4 _ 1 2 4 

2 y 4 2 



1-1/4 
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ig/, 2 + \z^ - \te 



1/2 



(61) 



where z < zi, i.e. z 2 h A /4 > c 2 z A /2. In using these expressions it has to be 
remembered that the moduli of the integrals have to be taken. To the right 
of the turning point at Z\ these solutions match on to 



1/wkb( 2 ) 



I o/i 2 _ l _ z 2 h A + I c 2 4 

2 y 4 2 



-1/4 



X cos 



l qh * _ 1^4 + I 2^4 

2 y 4 2 



1/2 



2/WKB ( Z ) 



1^2 _ 1^4 I 2^4 

2 y 4 2 



-1/4 



x sin 



ign 2 - ±*W + ~cV 



1/2 



**Ref. [22], p. 291, equations (21), (22) or Ref. ESI, Vol. I, Scc.6.2.4. 
ttRef. E3. P- 29L 



(62) 
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We now come to the algebra of evaluating the integrals in the above 
solutions. We begin with the exponential factors occurring in Eqs. floT]) . i.e. 



E ± 
exp \ ± 



dz 



2 y 4 2 



n l/2 



exp 



exp 



exp 



± 



8c 2 



d 



=F 



8c 2 



2cV 

2cV 
_ ^ 4_ 



1 - 



2cV ^ 1/2 



/i 1 

3/2 2:1 „ u 



2g 



=F 



/? 6 2 
8c 2 3 



2cV ^ 3/2 



Z? 1 



=F - 
2 



gZ^ 

8c 2 

«1 



Z 2 Zi 2 [1 
1 j( 2c 2 z 2 ■ 



2c 2 z 2 ] 



1/2- 



h 4 



[1 - ^]V2 



3 



[^- 2 2]1/2^2C 2 ; 



1/2' 



(63) 



Here the first part is the exponential factor contained in yA{z),y A (z) respec- 
tively (cf. Eq. pijl ). In the remaining factor we have (looking up Tables of 
Integrals) 

dz 



Aw ~ z2 V' 2 



V h 4 



/2c s 



2\ 1/2 



2c 2 

^4Xl/2 

2c s 



1/2 



2c 2 



1/2 



2c 2 



1/2 



Since we are interested in determining the proportionality of two solutions in 
their common domain of validity we require only the dominant z- dependence 
contained in this expression. We obtain this factor by expanding the expres- 
sion in powers of z/z\ (since in the integral z < Z\). Thus the above factor 
yields 



/2c 2 
U 4 



2\ 1/2 



hi 



z\2c 2 ) 



1/2 
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so that (cf. Eq. (gOJ) 

E± = exp 

_ ^±9/2 



h 6 2 r 

1/2- 



-C- 



2c^ 

T9/2 



2^2 ^ 3/2 



4 \ 1/2 



( 2 ^) 



,±9/2 



cxp 



±i 



_ i z 2^4 1 2 4 

4 2 



1/2' 



(64) 



Thus at the left end of the domain of validity of the WKB solutions we have 



(Mi) / \ 

yWKB\ Z ) 



\z 2 W 



1/4 



and y§£l(z 



E_ 



\z 2 h* 



1/4' 



(65) 



Comparing these solutions now with the solutions ()20a|) and (j20b|) . we see 
that in their common domain of validity 



where 



or 



Va(z) = PVvnaiiz), Va{z) 



and (3 



2h 2 



~h 2 ~ 


1/2 


L (h 2 ) - 


q/2 


_ 2 _ 




. (2c 2 y/ 2 _ 





' h 2 ' 


1/2 




-q/2 


2 _ 




. (2c 2 )V2_ 









4)9/2 



(66) 



(67a) 



(67b) 



(2c 2 )V2_ 

apart from factors [1+0(1/ h 2 )]. In these expressions we have chosen the signs 
of square roots of h A so that the conversion symmetry under replacements 
q — > —q, h 2 — > —h 2 is maintained. 

Returning to the even and odd solutions defined by Eqs. (J21)) we now 
have 



y±w 



[yA(z)±y A {z)} 



(6* 
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Now in the domain z — > oo we have 



Zl L 



2 H 4 2 



1/2 



CZ 2 

dz—= 
V2 



Ac 2 z 2 



cz 



CZ" 



h 6 



3V2 12c 2 ' 



(69) 



Inserting this into the solutions (}6*2~|) and these into ()68|) we can rewrite the 
even and odd solutions for 9ftz — > oo as (by separating cosine and sine into 
their exponential components) 



1 2 4 


-1/4 




2 CZ _ 




S + (±)exp 



C2" 



/i 6 



-^-(i) exp 



\3Vl _ 12^ 



(70) 



where 



S + (±) 

sl(±) 



P±-p)exp[i 



71 



PT~p] exp ( 



(71) 



Imposing the boundary condition that the even and odd solutions have the 
asymptotic behaviour given by Eq. (|59jh we see that we have to demand that 



S+(±)=0, i.e. Ip±-P = 0. 

2 i 



(72) 



Inserting expressions ()67a|) for P and p, this equation can be rewritten as 



'-h 2 r' 2 



2 /2 iP _ 2^(-l)U h^\ 
] 2P~ v^T W) 



6 \ 5/2 



(73a) 



Thus we impose this second boundary condition by making in Eq. (156)) the 
replacement 

(_ /i 2 ) 50/2 > 



)50-l 



2c 2 



(73b) 
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Inserting this into the latter equation we obtain (the factor "i" arising 
from the minus sign on the left of Eq. (|73b[l ) 

with g = 1, 3, 5 . . .. 

3.6 The complex eigenvalues 

We now return to the expansion of the eigenvalues, i.e. Eq. 



Expanding about q = go we obtain 

E(q,h 2 ) = £ (g ,/i 2 ) + (g-g ) 

= E (q ,h 2 ) + (q-q y— + -- - . (75) 



h 2 



2 

Clearly the expression for (g — go) has to be inserted here giving in the 
dominant approximation 

1,6 \ w/2 
2 i°h 2 {^) h 6 



^^^^wfe 8 6c2 - (76) 

The imaginary part of this expression agrees with the result of Bender and 
Wu (see formula (3.36) of Ref. [3]) for h = 1 and in their notation 

qo = 2K+l, — = 6. 

For comparison with the case of the double well below we note here that the 
ratio (3/(3 which we required for the derivation of the result is the ratio of 
matching coefficients. Hence the result does not involve a specific normali- 
sation of the WKB solutions. This is different in the case of the double well 
potential that we consider below. 
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4 The Double- Well Potential 



4.1 Defining the problem 

In dealing with the case of the symmetric double-well potential, we shall 
employ basically the same procedure as above. But there are significant 
differences. 

We consider the following equation 

^ + [E — V(z)}y(z) = (77) 
with double-well potential 

V{z) = v(z) = ~\z 2 h A + \c 2 z A for c 2 > 0, h A > 0. (78) 




Fig. 4 The double-well potential. 



The minima of V(z) on either side of the central maximum at z = are 
located at 

* = 4c < 79) 
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with 



V(z ± ) 



h 8 



2 5 c 2 ' 

h\ V {3 \z±) = ±6ch 2 , 
12c 2 , V {i) (z±) = 0, i > 5. 



In order to obtain a rough approximation of the eigenvalues we expand 
the potential about the minima at z± and obtain 



dz 2 



+ 



E - V(z ± ) --{z- z ± fh A + 0[(z - z ± ) 3 ] 



y = 0. 



We set 



-y/2h 2 

(thus we sometimes use /i 4 and sometimes h±) and 



E - V(z ± ) = l -q ± hl + A 



With the further substitution 



lj± = h±(z — z±] 



Eq. (jSTjl becomes 



where 



^9±(^±)2/(^±) = 



4 



(82) 



(83) 



(84) 



(85a) 



2^' (85b) 
By comparison of Eqs. ()85a|) . (J85bj) with the equation of parabolic cylinder 
functions u(z) = DJz), i.e. 



d 2 u{z) 
dz 2 



+ 



v H 



u(z) = 0, 



we conclude that in the dominant approximation q± is an odd integer, q = 
2n + l,n = 0,l,2,.... Inserting the expression (|83|) for E into Eq. ()77|) we 
obtain 



dz 2 



1 9 A 1 , , , 
2^ + 74-4^ 



(86) 
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where 



U{z) = vr [V{z)-V{z ± )], 



(87a) 



and near a minimum at z± 

U{z) = {z-z ± ) 2 + 0[{z-z ± f}. 



(87b) 

Our basic equation, Eq. ([8*BJl. is again seen to be invariant under a change of 
sign of z. Thus again, given one solution, we obtain another by replacing z 
by -z. 

4.2 Three pairs of solutions 

We define our first pair of solutions y(z) as solutions with the proportionality 

1 



y{z) oc exp 



U 1/2 {z)dz 



(88) 



Evaluating the exponential we can define these as the pair 
y A (z) = A(z)exp 
Va( z ) = A(z)exp 



v^\3 4c" 

i r c 



h 4 



(89) 



The equation for A(z) is given by the following equation with upper signs 
and the equation for A(z) by the following equation with lower signs: 



A"(z) T v^jcz 2 - ^A'(z) T V2czA(z) + {^-q ± h 2 ± + A 



= 0. (90) 



Since 



2+ 



2c' 



hi = V2h 2 , 



and selecting z + with g + = g, these equations can be rewritten as 



z 2 )A\z) + (qz + - z)A(z) = 



A"(*) + ^4(z) 



(4 - z 2 )A (z) - ( 9 z + + z)A(z) 



V2 
2c 



+ ^A{z) 



(91a) 
(91b) 
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To a first approximation for large h 2 = 2cz + we can neglect the right hand 
side. The dominant approximation to A is then the function A q given by the 
solution of the first order differential equation 

d 

dz 

We observe that a change of sign of z in this equation is equivalent to a change 
of sign of q, but the solution is a different one, i.e. A q (z) = A q (— z) = A^ q (z). 
Integration of Eq. (}9*2"j) yields the following expression 



V q A q {z) = 0, V q = (4 - z 2 ) — + {qz+ - z) . (92) 



A q (z) 



1 


z — z + 


\z 2 — z\ 


1/2 


z + z + 



q ' 2 \z-z+\^ 



\Z + Za 



(93) 



Looking at Eqs. (|91aj) . (|91b|) we observe that the solution y A (q, h 2 ; z) may 
be obtained from the solution y^(q, h 2 ; z) by either changing the sign of z 
throughout or — alternatively — the signs of both q and h 2 (and/or c), i.e. 

Va(<1, h2 'i z ) = VA(q, h 2 ; -z) = y A (-q, ~h 2 ; z). (94) 

Both solutions yA(z),y A (z) are associated with one and the same expansion 
for A and hence E. We leave the calculation of A to Appendix C. The result 
is given by Eq. ([(J.lj) . i.e. 

A = -c 2 (3g 2 + l)-^g(17g 2 + 19) + ---, 

jp ( h 2, h* 1 c 2 (3g 2 + l) 

E{q,h) = --— + — =qh 



2 5 c 2 y[T 2h A 
!C -g(17g 2 + 19) + of-^) 



The solutions y A (z),y A (z) derived above are valid in the domains away 
from the minima, 

z - z± I > O f 



hi 

We can define solutions which are even or odd about z = as 
y±{z) = ^[y A {q,h 2 ;z) ±y A (q,h 2 ;z)} 

= 7^[yA(q,h 2 ;z)±y A (q,h 2 ]-z)) 

= ^[y A (q,h 2 ;z)±y A (-q,-h 2 ;z)]. (96) 
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Considering only the leading approximations of the unnormalised solutions 
considered explicitly above, we have (since A q (0) = l/z + = A q (0),A' q (0) = 

-?/4) 

y + (q,h 2 ;0)='^, y_(q, h 2 ; 0) = 0, 
y' + (q, h 2 - 0) = 0, y'_(q, h 2 - 0) = A_ - (97) 
Our second pair of solutions, yB(z),y B (z) is obtained around a minimum 



of the potential. We see already from Eqs. (j85aj) and (|85bJ), that the so- 
lution there is of parabolic cylinder type. This means, in this case we use 
the Schrodinger equation with the potential V(z) expanded about z± as in 
Eq. (IHTj) . Inserting (jSHJ) and setting u± = h±(z — z±), the equation is — with 
differential operator T> q as defined by Eq. (|85b|) ) — 



V q (u ± )y(u±) = 



± 2 5 / i ch 3 u 3 + + c 2 ui - 2 A 



y(^±). 



Thus we can write a first solution 



/ 1 \ D l(a-l)( w ±( z )) 
y B {z) = B q [w ± (z)]+0{^y B q [w ± {z)] = ^_ 1 — , (98a 



[I( ? -l)]!2iM 



and another 



Vb{z) = Vb(-z) = B q [w±(z)}+o(J^\ 

= B q [w ± (-z)] + o(J^. (98b) 

Again the higher order terms along with the eigenvalue expansion in terms 
of q are obtained perturbatively. 

It is clear that correspondingly we have solutions yc(z),y c (z) with com- 
plex variables and C q {w) given by Eq. (J23J) with appropriate change of pa- 
rameters to those of the present case. Thus 

y c (z) = C q [w ± {-z)]+O f 



C q [w ± {-z)] = + , (99a) 
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y c (z) = y c (-z) = C q [w±{-z)\ = C q [w±(+z)] + of-^j . (99b) 

These are solutions again around a minimum and with yB{z),y c {z) P ro ~ 
viding a pair of decreasing asymptotic solutions there (or correspondingly 
y B (z), yc(z)). We draw attention to two additional points. Since w±(—z) = 
—h±(z + z±), we have w±(— z±) = —2h±z±, but w±(z±) = 0. Moreover, in 
view of the factor "i" in the argument of C q the solutions y^,y c have the 
same exponential behaviour near a minimum. 



4.3 Matching of solutions 



Next we consider the proportionality of solutions yA{z) and ys{z). Evaluat- 
ing the exponential factor contained in Ua{z) of Eq. (|89|) for z — > z±, we have 
(cf. (TO ) 



exp 



U xj2 {z)di 



exp 
exp 
exp 



-hi 



-hi 



z — z±)dz 



-z — zz± 



-hl(z-z ± f 



exp 



-h 2 z 2 



Allowing z to approach z + in the solution yA(z), we have 



i \-(a-i) ( w +\Ug-i) 

o l*-*+l 2W _ e R4 fi -fo+(«-«+) 3 = Kh+> . -^* 2 - 1 ™ 8 



|2zJ 5(9+1) 



(2z+) 



K9+1) 



-e 4 



e 4 



. (100) 



Recalling that around axgw± ~ 0, the dominant term in the power expansion 
of the parabolic cylinder function for large values of its argument is given by 

D v {w ± ) ~ wle- w ±'\ 

and comparing with Eq. ()98aj) we see that (considering only dominant con- 
tributions) in their common domain of validity 



yA{z) = -yB\z , a 

a 



(h + )^ q - l \2z + )^+ l ) 
2^-D^-i)]! 
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-±h 2 Z 2 



1 + 0| « 



. (101) 



Similarly we obtain in approaching z + \ 



Va{z) 



(2z + ) 



1(9-1) 



[z - Za 



(9+1) 



(102) 



and 



D 



Vc\ z ) 



))2i(" +1 ) 



2i(3+ 1 ) 



1)]! 



i-{(q+l)]\[ih + (z-z + p(«+V 
Therefore in their common domain of validity 



(103) 



a 



2i(" +1 ) 



h- 



Q 



(2z, 



,1(9- 



l) (-hl 



I X (9+1) [ 



!(?+!)]! 



1 + 



1 

hi 



(104) 

We have thus found three pairs of solutions: The two solutions of type A 
are valid in regions away from the minima, and are both in their parameter- 
dependence asymptotically decreasing there and permit us therefore to define 
the extensions of the solutions y± which are respectively even and odd about 
z = to the minima. The pair of solutions of type B is defined around 
argz = 0,7r and the solutions of type C around argz = ±7r/2. The next 
aspect to be considered is that of boundary conditions. We have to impose 
boundary conditions at the minima and at the origin. The solutions in terms 
of parabolic cylinder functions have a wide range of validity, even above the 
turning points, but it is clear that none of the above solutions can be used 
at the top of the central barrier. Thus it is unavoidable to appeal to (other 
and that means) WKB methods to apply the necessary boundary conditions 
at that point. The involvement of these WKB solutions leads to problems 
since, basically, they assume large quantum numbers. Various investigations, 
such as those of Refs. [T7] , [TH] , [IH] , [21] , therefore struggle to overcome this 
to a good approximation. We achieve the same goal (i.e. approximation) 
as such corrections here by demanding our basic perturbation solutions to 
be interconvertible on the basis of the parameter symmetries of the original 
equation. 



4.4 Boundary conditions at the minima 

(A) Formulation of the boundary conditions 
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The present case of the double-well potential differs from that of the 
simple harmonic oscillator potential in having two minima instead of one. 



Since it is more probable to find a particle in the region of a minimum than 
elsewhere, we naturally expect the wave function there to be similar to that 
of the harmonic oscillator, and this means at both minima. Thus the most 
basic solution would be even with maxima at z±, as indicated in Fig. 5. 
However, an even wave function can also pass through zero at these points, 
as indicated in Fig. 5. The odd wave function then exhibits a correspondingly 
opposite behaviour, as indicated there. At §lz — > ±00 we require the wave 
functions to vanish so that they are square integrable. We have therefore the 
following two sets of boundary conditions at the local minima of the double 
well potential: 






Fig. 5 Behaviour of fundamental wave functions. 




(105) 



and 




y+0±) 




(106) 
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We have 



y±{z±) 



[y A (z)±y A (z)] z _ 



-y B {z±) ± -y c {z±) 
a a 



and 



y'±{z±) 



-y' B {z±) ± = tfc( z ±) 
a a 



Hence the conditions ()105J) . ()106|) imply 



Vb(z±) = and VbM = « 

Vc( z ±) «' S7c(*±) a" 



(107) 



(10? 



(109) 



(B) Evaluation of the boundary conditions 



Inserting into the first of Eqs. ()109|) the dominant approximations we 
obtain (cf. also Eqs. (lETTaJ) and <IE7TE1> ) 



T (_l)|(5+1) 



2^[\{q-l)\\ 



-e 2 ' 



=F 



2 9/2 [|(g-l)]![l(g-3)]!sin{f(g + l)} 
1 J V2[i(g-l)]!sin{f(g + l)} 



(110) 



where we used first the reflection formula and then the duplication formula. 
Thus 



ain{j(g+l)} = T(-l)i^^ 



7T (/&)«/ 2 (2z+)« _i. 



(9-1)]! 



-e 2 



(111) 



Using formulae derived in Appendix B we can rewrite the second of 
Eqs. dnSD as 



^ = ^{I( 4 -S)}^cat{l( 4+ 1) 



a 

=F=- 

a 



(112) 
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Using Eq. (|11U|) this equation can be written 



-zcos<| 1) \ = =F 



Now 



sin<J^(g + l)| ~ sin|^(g + 1)| + ^(g-g )cos|^(g 
~ (-ljito+D^-^E for go = 3>7jll> 



(114) 



and 



cos|^(g + l)j ~ cos|^(g + 1)| - ^(g-go)sin|^(g + 1) 

~ _(_1)*(«>-I)(g_g )?[ for g = l j5)9) ... 

= -(g-g )^(-l)^ 0+1) (-l) 1/2 - (H5) 

Thus altogether we obtain from the boundary conditions at z + (and corre- 
spondingly from those at zJ) 

In Appendix E (after calculation of turning points) we rederive this relation 
using the WKB solutions from above the turning points matched (linearly) 
to their counterparts below the turning points and then evaluated at the 
minimum. 

In summary: We needed the solutions of type A for the definition of even 
and odd solutions. Since the type A solutions are not valid at the minima, 
we matched them to the solutions of types B and C which are valid there 
and hence permit the imposition of boundary conditions at the minima. 



4.5 Boundary conditions at the origin 

(A) Formulation of the boundary conditions 
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V(z) 




Since our even and odd solutions are defined to be even or odd with 
respect to the origin, we must also demand this behaviour here along with a 
nonvanishing Wronskian. 

Fig. 6 Turning points zq,z±. 
Hence we have to impose at z = the conditions 

j/_(0) = 0, y'_(0)^0, 

y+(0) ± 0, j/ + (0) = 0. (117) 

Thus we require the extension of our solutions to the region around the local 
maximum at z — 0. We do this with the help of WKB solutions. We deduce 
from Eq. (|97|) that the two turning points at z and Z\ to the right of the 
origin are given by 



\qhl + ^ _ 1 h l U{z) = o, (118) 

i.e. 



\ a h 2 + -z 2 h A - -c 2 z 4 - — = o(— 
2 g+ + 4 2 2V U Ul 
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Using z + = h 2 /2c, one finds that 

I- 

\4c 2 



2c 



(2g 



2\l/4 



/l 



and 



4c 2 



+ 



qh\ 



+ 



A 
hi 



2 1 / 2 cg 

1/2 



A 



1/2 



2c 



(2g 



2U/4 



(119) 



(120) 



The height of the potential at these points is 



V{z)\ 



20, Zl 



h s qh 
+ 



2 5 c 



2V2- 



Thus for large values of h 2 the turning points are very close to minima of the 
potential for nonasymptotically large values of q. 

(B) Evaluation of the boundary conditions 

We now proceed to evaluate the boundary conditions (|117|) . In the domain 
< z < Zq, i.e. to the left of z where V > E, the dominant terms of the 
WKB solutions are* 



Vwkb( z ) 



- \qh\ + \h%U{z) 



1-1/4 



x exp 



dz 



qh\ + \h%U{z) 



1/2 



'121) 



and 



- \qh\ + \h%U{z) 



1-1/4 



x exp 



■-0 



dz 



- \qh\ + \h%U{z) 



1/2 



(122) 



In order to be able to extend the even and odd solutions to z — 0, we have to 
match yA.(z),y A (z) to ^wKm^) an< ^ ^wkb( z )- We therefore have to consider 



The superscript (I, zq) means "to the left of zq" 
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the exponential factors occurring in (|121|) and (J122J1 and consider both types 
of solutions in a domain approaching but not reaching the minimum of the 
potential at z + . Thus we consider in the domain \z — z+\ > (2q/h 2 + ) l l 2 : 



h 



-0 



dz 



4 k * 



\qh\ + \h%U{z) 



1/2 



dz 



( ^2 



20 — 2+ 



2q 



nl/2 



dz 7 z'- 



1/2 



"IT 



1/2 



5 1 
_ /d ln 



,/2 



2q 
hi 



1/2 



20-2+ 



2 — Z+ 



(123) 



Evaluating this we have (with z — z + = — {2q 2 ) l / 4 /h = — {2q/h 2 + ) 1 / 2 ) 



1.2 W 2 1 



±/l 



/2 



h\z> 2 

1/2 



z'=-(2g/hi ) 1 / 2 



1/2 



1 1 , 



z — z. 
h 2 



2g 
hi 



1/2 



2? 



~/^*-* + ) 2 + -gln|2(z-z + )|. 



(124) 



In identifying the WKB exponentials we recall that ^wkb is exponentially 
increasing and ?/wkb exponentially decreasing. Hence we have for ?/wkb the 
exponential factor 



cxp 



20 



dz 



\qh\ + \h%U{z) 



1/2 



-z + \^eh(\hl 



g/4 



-9/4 



-\{z-z + )*h\ 



(125a) 
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We observe here that the WKB solution involves unavoidably the quantum- 
number-dependent factors exp(g/4) and (g/4) 9 / 4 which do not appear as 
such in the perturbation solutions. The only way to relate these solutions 
is with the help of the Stirling formula which converts the product or ratio 
of such factors into factorials such as those inserted from the beginning into 
the unperturbed solutions (|98a|) and (|99a|) . However, Stirling's formula is 
the dominant term of the asymptotic expansion of a factorial or gamma 
function and thus assumes the argument (oc q ~ In + 1) to be large (it is 
known, of course, that the Stirling approximation is amazingly good even 
for small values of the argument). Thus, using the Stirling formula z\ ~ 
e -(*+i)( z + 1)^+2^/2^, we can rewrite the exponential as 

I l 1 V2n 

exp ( - J dz - -qh 2 + + -h\U {z) 

\y — y 19/2 /I \ ?/ 4 

- (2-) 1/2 | 7 ^ ± | jI y(^l) e-K»+r*. (125b) 

Here q/4 was assumed to be large but we write [j(q — 1)]! since this is the 
factor appearing in the solution (|98a[) . We see therefore, since there is no way 
to obtain an exact leading order approximation with Stirling's formula for 
small values of q, the results necessarily require adjustment or normalisation 
there in the g-dependence. This is the aspect investigated by Furry 
Since correspondingly 

1 1 1" 1/4 2 1 / 2 

-qhl + Z h + U{z)\ ^ {z _ z+y/ 2 {hi y^ 

we obtain 

VWKBW ~ ^4)1/4 [l( g _l)]| 

2q- 
M 

This expression is valid to the left of the turning point at zq above the 
minimum at z + . 

In a corresponding manner — i.e. using Stirling's formula (and not the 
inversion relation) — we have 



1/2 



for \z-z+\ > [ -rt J . (126) 



0F(/4)V4| z _ z+ |!( 3 +i) V2 
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where was written as [\{q — 3)]! for q large. Comparing for z — > z + 
the WKB solutions ljT2K|) . (|T2T|) with the type-A solutions (HUDJ), (jTH^I we 
obtain in leading order (i.e. multiplied by (1 +0(l/ti+))) the proportionality 
constants 7, 7 of the matching relations 

= = 1Va{z): (128) 



i.e. 



and 



7 



20F A 



(^ 4 + ) 1/4 [J(?-l)]! 



_ \w-mfi 

7 



/1 \ 9/4 

[^h%\ {2z + )^ q+1) e- 1 * h + z2 + (129a) 
1 \ ~ 9/4 

^] (2^+)-3(*- 1) ei fc +*+. (129b) 



0F(^ + )V4 

Using again the duplication formula''' the ratio of these constants becomes 

7 .^yw. K< 



7 [Ks - 

Since the factorials [j(q — 1)]!, — 3)]! are really correct replacements of 
[|(?)]! only for q large, this result is somewhat imprecise. However, it is our 
philosophy here that the factorials with factors occurring in the perturbation 
solutions are the more natural and hence correct expressions, as the results 
also seem to support. The relation ()130|) is used in Appendix E for the 
calculation of the tunneling deviation q — q by using the usual (i.e. linearly 
matched) WKB solutions, and is shown to reproduce correctly the result 
()116|) which was obtained with our perturbation solutions from the boundary 
conditions at a minimum. 

Returning to the even and odd solutions f)9fij) we have 

V±W = \[va{z) ±y A {z)\ = ^y^l(z) ± ^y^l(z). (131) 

Applying the boundary conditions ()117j) we obtain 

yS(o) 7 j#S'(o) 7 



(132) 



^In the present case as 



~ - 3)]! = v^2-^9-D[I ( g _ !)],. 



42 



Thus we have to consider the behaviour of the integrals occurring in the 
solutions (QUI), (HID near z = 0. We have 



-5(1 



l -qh\ + ^ 4 + [^) 



1/2 



2(] 



dz 



1/ 



_ _ 1 2, 4 , 1 2 4 , 

2 q + 4 + 2 + 25 C 2 



1 1/2 



1 fh 4 



2 \2 2 c 



cz 



2 1 + 



"2 

— = I dz 

V2 



h 4 



1/2 



fh A 



__ >/ g_i_^)[_ + v ^_±_ Z 



Setting 



/) 4 7) 

= + ~ 4, a 2 



6 =4?-^c 



\4c 2 



1/2 



(133) 



we can rewrite the integral as 



-= / dzyj{a 2 -z 2 )(b 2 -z 2 ). 
v2 J z 



The integral appearing here is an elliptic integral which can be looked up in 
Ref. [2*3 The elliptic modulus k (with k' 2 = 1 — A; 2 ) and an expression u 
appearing in the integral are defined by 



k =~ 

a z 



1-u 



l + u h 3 



u = TXv 7 ^ = Gy 2g, G 



8y/2, 



c 



so that 



(134a) 
(134b) 



The integral I^iz) evaluated at z = is then given in Ref. [2H] by 



ca 



[(1 + k 2 )E{k) - k' 2 K{k)} 



nVzc 2 

2 



y/T+^[E(k) -uK(k)}, 



_3G 2 



VT+u[E(k) -uK(k)], 



tSee Ref. formulae 220.05, p. 60 and 361.19, p. 213. 



(135a) 
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where K(k) and E(k) are the complete elliptic integrals of the first and 
second kinds respectively. Here we are interested in the behaviour of the 
integral in the domain of large h 6 /c 2 which implies small G 2 . The nontrivial 
expansions are derived in Appendix F, where the result is shown to be (with 
g~g = 2n + l,n = 0,l,2,...) 



h(0) 



T A ± 2V2 12c 2 



h 6 q 
^ - In 



2h 3 



"3G 2 



2 1 /4 cg V2 



G 



i n + - In - i- n+- In n+- 



(135b) 



in agreement with Ref. JD]^ Since the integral is to be positive (as required 
in the WKB solutions) we have to take the upper signs. It then follows that 
(again in each case in leading order) 



2/wkb(^)I«=o 



(136) 



2 5 c 2 ' V4' 

where h 8 /2 5 c 2 = —V(z±). Using Stirling's formula we can write this 



^WK°b( Z )Iz=0 

Correspondingly we find 



2/wkb(^)U=o 



23/2. 



(2&) 1/4 [|(?-1)]! 



g \2\f2c l 



(137) 




\iiQ-m 



and 



and 



dz 



Vwkb( z j 



2 = 



h s y/4 (27r) i/ 2( ^_^ )g /4 



g 12v / 2c 



2 5 c 2 



e 12 



V2^ 



d_ 

dz 



2=0 



2 5 c 2 



1/4 



[|(9-3)]! 
(2tt)V 2 (^^) 9 /4 



g 12v^c 2 _ 



;i38) 



(139) 



(140) 



^The expansion of I2 used in Ref. 1 2 1J I misses an n-dependent power of 2 in the result. 
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With these expressions we obtain now from Eqs. ()132|) — on using once again 
the duplication formula in the same form as above — 



^SFlM e ^ = ? (141) 

Comparing this result with that of Eq. (J13Uj) . we can therefore impose the 
boundary conditions at z = by making in Eq. (J116J) the replacement 

{h\yl\2z^e~-^l ^l^j^j e~W. (142) 

A corresponding relation holds for h + and z + replaced by /i_ and z— In- 
serting the expressions for h + and z + in terms of h, we see that the relation 
is really an identity (the pre-exponential factors on the left and on the right 
being equal) with the exponential on the left being an approximation of the 
exponential on the right (which contains the full action of the instanton). 
This is our second condition in the present case of the double well potential. 



4.6 Eigenvalues and level splitting 

We now insert the relation ()142j) into Eq. (jllfij) with go = 2n + 1, n 
0, 1,2,.. ., and obtain 



{q - *> * W ^w^;-!)]! 6 ^- * = 3 > 5 - • • • ■ (i43 » 

We obtain the energy E(q, h 2 ) and hence the splitting of asymptotically de- 
generate energy levels with the help of Eq. (|1()6J) . Expanding this around an 
odd integer q we have 

/dE 
V 9q 



qo 

h 2 



E{q,h 2 ) ~ E (q ,h 2 ) + (q-q )i 

~ E (q ,h 2 ) + (q-q )^=. (144) 
Inserting here the result ()143|) we obtain for E(q , h 2 ) the split expressions 

E ± (q , h 2 ) ~ E (q , h 2 )T ^^rf^Z l)]i e ~^' 96 = 1,3,5,..., (145) 
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where E (q , h 2 ) is given by Eq. (jC.l|) . i.e.' 

h 8 1 , 2 
'2*c 2 + ^2 2h* ' ~\h™ 



r, / ,2, h 8 1, 2 c 2 (3q 2 + 1) J 1 



Thus 



AE(q 0: h 2 ) = E4q Q ,h 2 )-E + (q ,h 2 ) 

2 *o+2 h 2 / ^6 \ 90/2 



e 2172^2 I m ass m = - 
y/*2*>/4[±(q -l)]\\2c*J V 2 

2 (2n+9)/4 ft (M g-ipfo. (146) 



Combining Eqs. (fTlij) and ()E.14|) with (fT32Jl . the level splitting, i.e. the 
difference between the eigenenergies of even and odd states with (here) go = 
2n + 1, n = 0, 1, 2, . . ., for finite /i 2 , can be given by 

(147) 



l( d_E \ yjj^(0) 



In Ref. the WKB level splitting is effectively (i.e. in the WKB restricted 
sense) defined by 

A WKB _ 1 ^E 2/wKb(0) _ 1 a EV„ /,2\ (TAQ\ 

" "^(" + I)5S(0)"2 A£( *' M - (148) 

Here dE/dn corresponds to the usual oscillator frequency. The result (|146|) 
is in Ref. ^H] described as the "modified well and barrier 11 result A^ WB , the 
pure WKB result of Ref. [TH] (i.e. that without the use of Stirling's formula 
and so left in terms of e and n n ) being this divided by the Furry factor 



fn 



2vr \n+ \ 



2 



^In Ref. )18l the "usual WKB approximation" of this expression is — with replacements h 4 /4 
fc, c 2 /2 <-> A — given as 

i.e. the following expression for large q ~ 2n + 1, which evidently supplies some correction terms. 
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which is unity for n — > oo with Stirling's formula." Of course, these deriva- 
tions do not exploit the symmetry of the original equation under the inter- 
changes q <-» —q, h 2 «-> —h 2 , as we do here. We see therefore, that if this 
symmetry is taken into account from the very beginning (with the appro- 
priate use of the Stirling formula) the Furry factor corrected result follows 
automatically. The Furry factor represents effectively a correction factor to 
the normalisation constants of WKB wave functions (which are normally 
l/y/2n and independent of n as explained in Ref. [2T]) to yield an improve- 
ment of WKB results for small values of n, as is also explained in Ref. [TH].** 
Had we taken the mass mo of the particle in the symmetric double-well 
potential equal to 1 (instead of 1/2), we would have obtained the result with 
E, h 4 and c 2 replaced by 2E, 2h 4 and 2c 2 respectively (see Eq. (JjJJ)). Then 



E(q ,h 2 ) = E (q ,h 2 )T 



V^2W4[I (gb _ 1)] 



_ h 
■e e^ 2 



(149) 



and AE becomes 



A 1 E(q ,h 2 )~2 



2h 2 



7r2W4[I( go _i)]!^2V2 c 2 



L6 \ 9o/2 

e~ h6/6c \ (m 



with 



E {qo,h 2 



h 8 



\m =l 



2 4 c 2 



q h" 



If in addition the potential is written in the form 

A 



V(z) 



, 2 -^V, A>0, 



(150) 



(151) 



II Ref. HI], Eq. (66). 

"Comparing the present work with that of Ref. [181 . we identity the parameter k there with ft 4 /4 here, 



and A there with c 2 /2 here. Then the splitting Aj^ WB of Eq. (24) there is 



^"2( 10 " + 13 )/ 4 /fc3/2\ n + l/2 



A 



cxp 



3 A 



h 2 



2 (2n+5)/4 



/ h 6y l+ i/2 


h 6 1 


{*) CXP 


6a/2c 2 . 



agreement with AE(q , h 2 )/2 of Eq. HHH . 
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a form frequently used in field theoretic applications, so that by comparison 
with Eq. dZHJ) h A = 2/i 2 , c 2 = A/2, the level splitting is 



This result agrees with the ground state (n = 0) result of Ref. P using 
the path-integral method for the evaluation of pseudoparticle (instanton) 
contributions.* Equation ()150j) agrees also with the result of Ref. ^Hj for 
arbitrary levels obtained with the use of periodic instantons^ and with the 
results of Ref. [Hj using multi-instanton methods. * 

5 Concluding Remarks 

In the above we have presented a fairly complete treatment of the large- 
h 2 case of the quartic anharmonic oscillator, carried out along the lines of 
the corresponding calculations for the cosine potential and thus of the well- 
established Mathieu equation. In principle one could also consider the case 
of small values of h 2 and obtain convergent instead of asymptotic expan- 
sions; however, these are presumably not of much interest in physics. We 
considered only the symmetric two-minima potential. The asymmetric case 
can presumably also be dealt with in a similar way since various references 
point out that the asymmetric double well case can be transformed into a 
symmetric one (see e.g. Refs. pp. [TB]). 

Every now and then literature appears which purports to overcome the 
allegedly ill-natured "divergent perturbation series" of the anharmonic oscil- 
lator problem, and reasons and even numerical studies are presented to sup- 
port this claim. [3U] It is clear — as also demonstrated by the work of Bender 
and Wu [I] , [2] — that the expansions considered above are asymptotic. Ta- 
bles of properties of Special Functions are filled with such expansions derived 

*See Ref. 0, form ula ( 4.11). The definition as A^KB is usedi so that thc rcsu i t f Eq. differs 
by the factor 2 in Eq. fTISI . 

^To help the comparison note that in Ref. |10l the potential is written as 



for mo = 1. The comparison with Eq. 11491 therefore implies the correspondence r\ 2 <-> A/2,m 2 <-> /J, 2 /2, 
and g 2 is given as g 2 = r/ 2 /m 3 . 

*Ref. [TU, first paper, Eqs. (2.34) and (E.15). 




(152) 
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from differential equations for all the well-known and less well-known Special 
Functions. There is no reason to view the anharmonic oscillator solutions 
differently. In fact, in principle the Schrodinger equation of the quartic oscil- 
lator is an equation akin to equations like the Mathieu or modified Mathieu 
equations which lie outside the range of equations of hypergeometric type. 
The immense amount of literature meanwhile accumulated for instance in the 
case of the Mathieu equation can indicate what else can be achieved along 
parallel lines in the case of special types of Schrodinger equations, like that 
for anharmonic oscillators. Conversely, new methods discovered for dealing 
with the quartic oscillator could similarly be applied to the periodic potential 
and tested there. 

The ground state splitting of the symmetric double well potential has been 
considered in a countless number of investigations. A reasonable, though 
incomplete list of references in this direction has been given by Garg be- 
ginning with the well-known though nonexplicit (and hence not really useful) 
ground state formula in the book of Landau and Lifshitz [32] ■ Very illuminat- 
ing discussions of double-wells and periodic potentials, mostly in connection 
with instantons, can be found in an article by Coleman (23] - The wide pub- 
licity given to the work of Bender and Wu j2] , [2j made pure mathematicians 
aware of the subject; as some relevant references with their approach we cite 
Refs. HH, [SH], [SE!,[SI!, [SB], [3H]. The double- well potential, in both sym- 
metric and asymmetric form, has also been the subject of numerous numerical 
studies. As references in this direction, though not exclusively numerical, we 
cite papers of the Uppsala group jJQ]. Perturbation theoretical aspects have 
also been employed in Ref. j3T]. Wave functions of symmetric and asym- 
metric double- well potentials have been considered in Ref. |12]> i n which it 
is demonstrated that actual physical tunneling takes place only into those 
states which have significant overlap with the false vacuum eigenfunction. 
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Appendices 



A Calculation of eigenvalues with the solu- 
tions of type A (inverted double well) 



Here we use the solutions of type A, i.e. (|20aj) . (|20b|) . to obtain in leading 
order the eigenvalues E, i.e. 



'-Aq 2 



o 



(A.r 



We rewrite the upper of Eqs. (|13|) in the following form with power ex- 
pansion of the square root quantities and division by h 2 : 



-zA'(z) + -(q-l)A(z) 



+ £ ( 

i=l x 



2, v 2\ 1 



2c 2 z 



\aizA\z) 



where the expansion coefficients are given by 



(A.2) 



a = l,«i = — -, a 2 



Po = 1, Pi 



a 3 



16 : 



«4 



128 



fa 



03 



7_ 
16' 



Using Eqs. (fT7j) and (fTHjl we obtain 



and 



i (g-i) 

2 z 



i(g-l)(g- 3)4^0*). 



(A.3) 
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The lowest order solution = A q (z) therefore leaves uncompensated on 
the right hand side of the equation for A the terms amounting to 

^ 0) = -^- 1 ^- 3 ^-*-w Aq ~w iq+2)Aq+4 



[(q - lfa + Pi]A q . 



i=3 

Clearly one now uses the relation 

z 2i A q (z) = A q+M (z). 

In this way Rq is expressed as a linear combination of functions A q+ n(z). 
As always in the procedure, one observes that 

V q :=-z^ + ^(q-l), V q+4i = V q + 2i, V q ^^ = fiA q+4l . (A.4) 

Thus a term [iA q+4 i in R q can be taken care of by adding to A^ the 
contribution / ^+ 4 ' except, of course, when i — 0. In this way we obtain the 
next order contribution A^ to A^ and the coefficient of terms with i — 0, 
i.e. those in A q (z), give an equation from which A is determined. Hence in 
the present case 

Am = __L fa _ 1)(? _3)^-^( ? + 2)^i- 

In its turn A^ leaves uncompensated terms amounting to 

= ~(Q 3)^{ (- ^) (?-5)(?-7)V. 

c 2 A f2c 2 \ 2 a 

(Q ~ 2 ) A i - - TT 77^+4 



A 



2/i 6 



A 9+4 + • • • > + 
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The sum of terms with A q in R q °\ R^q \ ■ ■ ■ must then be set equal to zero. 
Hence to the order we are calculating here 



(4/z 2 )(4/i 4 , 



[(q -l)(q- 2)(q - 3) - (q + l)(q + 2)(q + 3)] 



(A.6) 



i.e. 



It follows that 



A 



3 2/2 ^ « ' 1 



The same result is obtained in connection with the solution of type B, as 
explained in the text. 

B Evaluation of ^(0), y c {0), y' B (0), y' c (0) 

We show that — with w = hz — the leading terms of the quantities listed 
are given by 



B q (0) 

CM 



d . . 



d 

dw 



C q (w) 



i [Hq-m 



.IX 



sin{-(g + l)}, 



2 7T 

sm{- g + 3 }, 

71 4 



>/^ sin {^(?- 3 )}- 



(B.la) 



In fact, with the reflection formula (—z)\(z — 1)! = ir/ sin-zr^, one finds that 



\--(a + 111! = - 

[ ^ )V [|(g-3)]!sin{f(g + l)} 



and hence £^ = 1. (B.lb) 
C,(0) 
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From the literature, e.g. Ref. 



D 



(9-1) 



(0) 



we obtain 



-1(9 + 3)]!- 



(B.2) 



Thus with the help of the reflection formula cited above: 



B g (0) 



%9-D 



(0) 



[\(q - I)]l2z(«-V [\{q -l)]![-|(g + 1)]! 
1 [|(g-3)]! 



i[^-D]' Sm{ 4 (g+1)} 



(B.3) 



and 



2^ 



2 7T 

-sm{-(g + 3)}, 



(B.4) 



Expressions for C q (0), [C q (w)] w=0 follow with the help of the "circuit rela- 
tion" of parabolic cylinder functions given in the literature as 

%9-i)M = e-^^D^i-w) 



+ 



From this relation we obtain 



D 



,(0) 



D |(9-D 



(0) 



2 -(9+DV^ V2[i(g-l)]!cos{f(g-l)}- 
Inserting from ()B.2|) . using the above reflection formula and the duplication 
formula ■ s /tt(2z)\ = 2 2z z\(z — 1/2)!, we obtain 



D 



(9+1) 



(0) 



2i(ff+D [ _l(g_l)]| 



From this we derive 



5(9+1) 



(0) 



7T 



2^(9+l) [ _l (g+ l) ]! [_l(g+l)]![l(g_l)]! 



[3(5-3)]! . 7T 

■sm{-(g + 1)}. 



V^[i(9-1)]! 



(B.5) 
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Similarly we obtain 



and 



[^HU= 2i( ^_ (B.6) 



[ glWjJ ° [_I( g +l)]!2-i(^) [-K« + l)] ! [j(9-3)] ! 

= -zy / |sin{^(g+l)} = z v /|sin{^(g-3)}. (B.7) 

C Calculation of eigenvalues with solutions 
of type A (double well) 

We show in conjunction with the derivation of solution that the leading 
terms of A and hence E are given by 

A = -cW + 1) - ^|r<i(17« 2 + 19) + ■ ■ ■ , 

vy ; 2 5 c 2 2/i 4 

!C -g(17g 2 + 19) + of-LY (C.l) 



The structure of the solution (jHSJ) for A q (z) is very similar to that of the 
corresponding solution in considerations of other potentials, such as periodic 
potentials. Thus it is natural to explore analogous steps. The first such step 
would be to reexpress the right hand side of Eq. ()91a|) with A replaced by 
A q as a linear combination of terms A q+2 i, % = 0, ±1, ±2, . . .. First, however, 
we reexpress A" in terms of functions of z multiplied by A g . We know the 
first derivative of A q from Eq. (|9*2*|) . i.e. 

A' q (z) = ^^A q (z). (C.2) 

Z Z i 
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Differentiation yields 



A q (z) 



(z 2 - 4) 2 



(z - qz + ) 2 - (z 2 - zV) + 2z(z - qz + ) 



Ag{z) [2z 2 - Azz + q + z 2 + (q 2 + 1)]. 



(z 2 - z 2 ) 2 



(C.3) 



We wish to rewrite this expression as a sum 

coefficient iA q+2 i (z). 

i 

We also note at this stage the derivative of the entire solution Ua(z) 
into account only the dominant contribution: 



V2 



z — z. 



[z 2 - z 2 



A q (z) exp 



LS: 



■z z 

V^l3 4c 



We observe some properties of the function A g (z) given by Eq. (|93|) : 



A q +2iA q+ 2j 



A 



q+2i+2jA q , 



A 



q+2 



A n 



z — z + 
z + z + 



A n 



A q+2 + A g _ 2 _z 2 + z 2 + 



A q +2 — Ag_ 2 



A n 



' z 2 — zV 



A n 



A q -2 

zz + 



(4zz + ) 5 



A, 



r* 2 - ^i) 2 



A 9+4 - 2A q + A 



9-4) 



taking 
(C.4) 

(C.5) 

(C.6) 
(C.7) 



From these we obtain, for instance by componendo et dividendo, 

z 

z+ 



A q+2 + A q 



A q+2 - A q 
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Aq+2\ f A q + 



A n 



A n 
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1 + 2- 
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A n 



+ 2- 
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Similarly we obtain 



z+ 
z 



A q+ 2 - Aq 
Aq+2 + A q 
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A n 



+ 2- 
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(C.8) 
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(C.9a) 
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and from this 



= 1-4- 
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q+2 
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g+4 



A„ 
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12- 



.4 



g+6 



A n 



+ 16 : 



.4 
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A, 



(C.9b) 



Hence with Eq. ffTTTj) 



[A g+4 - 2A q + A 9 _ 4 ] 



= T + {^) 2[Aq+4 ~ 2Aq + Aq - 4 
= ^(4^) [ V4 - 2 V2 + + 2A q+2 

Finally with Eqs. (f(T7jl and ()C.9bjl we obtain 



1.2^ + 2^ 

Aq Aq 



Aq+A + 



(CIO) 



A, 



z 2 A 



(z 2 - 4) 2 (z 2 - Z 2 ) 2 Z 2 

1 \ 2 1 
— J [A q+4 -2Aq + Aq_ A ]^ 



1-4- 



.4 



g+2 



A, 



A, 



- 12- 



A 



A n 



= (^j [Aq_ 4 -4Aq. 2 +6Aq-4Aq +2 +A q+A + ---}. (C.ll) 

Inserting expressions (jC.7|) . ()C40|) and ()C.11|) into Eq. ()C.3|) . we obtain* 

A" q = (JA i(q-l)(q-3)Aq_ 4 -4( q -l) 2 Aq„ 2 + 2(3q 2 + l)Aq 

-4(g + l) 2 A q+2 + ( q + l)( q + 3)A q+4 ■ ■ ■]. (C.12) 



Hence the first approximation of A, 



'/■ 



*The reader may observe the similarity with the corresponding coefficients in the simpler case of the 
cosine potential, cf. Ref. [5]. Infinite series like here for solutions of type A arise in the more complicated 
cases, e.g. in the corresponding treatment of the elliptic potential (Lame equation); cf. Ref. 1271 . 
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leaves uncompensated on the right hand side of Eq. (j91a|) the contribution 



Rf 



2c [ A « + hi Aq . 

7^[{q, q ~ 4)A_ 4 + (q, q- 2)A q - 2 
+ (q,q + 4)A q+A + (q,q + 6)A q+6 4 



J,q)A q + (q,q + 2)A q+2 
}, (C13) 



where the lowest coefficients have been determined above as 

(q,q T A) = (q T l)(q T 3), (q, g =F 2) = — 4(g =F l) 2 , 

(q,q) = 2(3g 2 + l) + 164^ = 2(3g 2 + l) + 4A. (C.14) 

It is now clear how the calculation of higher order contributions proceeds 
in our standard way. In particular the dominant approximation of A is 
obtained by setting (q, q) = 0, i.e. 



A 

E(q,h 2 ) 



-c 2 (3q 2 + l] 



1/ 



2 5 c 2 v/2 



c 2 (3g 2 + i; 
2h A 



Since 
we have 



DqAq — 0, T> q+ 2iA q+ 2i 



and V q+2 i = V q + 2iz^ 



A 



q+2i 



-2%Za 



A 



q+2i, 



except, of course, for i = 0. The first approximation A^ = A q leaves 
uncompensated on the right hand side of Eq. ()91a)l the contribution R q °\ 
Terms /jA q+ 2i in this may therefore be eliminated by adding to A^ ' the 
contribution A^ 1 ' given by 



A 0) 



V2c 



(q, q - 4) 



A,-4 + 



(q,q-2] 



H ; ^ 9 +2 t 



8z + y 4^ 
(9,9 + 4) 



.4 



5-2 



-4z, 



V4 + 



(C.15) 
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The sum A = + A^ then represents a solution to that order provided 
the sum of terms in A q in R q ^ and Rq is set equal to zero, where Rq is the 
sum of terms left uncompensated by A^ l \ i.e. 



Ri l) 



V2c \ 
2 3 h*J 



(<?, 9-4) d(0 ) 



R q-A + 



> g ~ 2 ) ff (D) , 



(C.16) 



This coefficient of A g set equal to zero yields to that order the following 
equation 







^2c 



2 3 /i 4 



(q,q-A)(q-A,q) 



(q,q-2)(q-2,q) , {q, q + 2)(q + 2, g) 



4z + 

(g,g + 4)(g + 4,g) 
-8z+ 



-4z H 
+ ••• 



which reduces to 



= 2(3g 2 



V2& 
2h e 



g(17g 2 + 19), 



(C.17) 



thus yielding the next approximation of A. 



D Calculation of eigenvalues with solutions 
of type B (double well) 

Equation (|C.1J) or (|C.17J) can also be obtained in conjunction with the so- 
lutions of types B or C. The initial step is to use the recurrence relation 
of parabolic cylinder functions D u (w) to reexpress the right hand side of 
Eq. ()85a|) as a linear combination of functions B q (w) defined by Eq. ()98aj) . 
Thus with the recurrence relation 

wD u (w) = D v+1 (w) + vD v - X {w) (D.l) 
and the expression (|98a| . i.e. 

Di, ^(w) 

BJw) = 2(q } \ (D.2 
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we obtain 

The procedure is then similar to that in the case of the solution of type 

A. 



E Recalculation of tunneling deviation using 
WKB solutions 



Here we determine the tunneling deviation q — qo of q from an odd integer 
go, above obtained as Eq. (jllfij) . by now using the periodic WKB solutions 
below the turning points. 

The turning points at Zq and z\ on either side of the minimum at z + are 
given by Eqs. ffTTTl and (THZU1) . We start from Eq. (THTTl) . i.e. 



1 



y±{z) 



[y A (z) ± y A (z)] 



1 (l,zo) I \ I ^ _(£,Z()) / \ 



(e.i; 



Different from above we now continue the solutions (in the sense of linearly 
matched WKB solutions) across the turning point at Zq in the direction of 
the minimum of the potential at z + . Then ((r, Zq) meaning to the right of Zq 
and note the asymmetric factor of 2)t 



y±K z ) 



1 ,2 1,4 



qh 2 + - -h%U{z) 



-1/4 







1 -2 sin 


fdz( 


1 7 


L J z \ 



X 



± = cos 

7 



\qh\ - 



1/2 



\qh\ - \h%U{z) 



1/2 



(E.2) 



See Ref. EU, p. 291, Eqs. (21), (22) or Ref. |2J0, Vol. I, Sec. 6.2.4. 
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We also note that 
d 

Tz y±{z) r ~ 



\qh\ - \h%U{z) 



1/4 



ri 




< —2 cos 


/ dzt 


1 7 


- J 20 * 



f dz(±qh%- \h%U{z)X'\\ 



=F = sm 

7 



20 



1/2 



(E.3) 



We now apply the boundary conditions ()1U5|) and ()1U6|) at the minimum z + 
and obtain the conditions: 
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— sm 
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20 
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±— cos 


7* 
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- •'ZO 



1/2 
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(E.4) 



and 



— cos 

7 

1 

=F = sin 

7 



2+ 



20 



( -g^ - -hXU(z] 



1/2 



2+ 



20 



( -qh 2 + - -h 4 + U(z) 
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Hence 



and 



tan 



cot 



2(1 



1/2 



=F 



_7_ 

27 



Z+ dz(\qhl-\h%U{z] 



1/2 



7T 

+ 4 



±2. 

27 



(E.6) 



(E.7) 



In the present considerations we approach the minimum of the potential at 
z + by coming from the left, i.e. from z . We could, of course, approach the 
minimum also from the right, i.e. from z\. Then at any point z G (zq, Z\) we 
expect* 

(E.8) 



I (r,2o)/ \ 



I?/wkb( z ) 



i_(r,zo) / \ 

I^wkbI^J 



I?/wkb( z ) 



*See e.g. Ref. H3, Sec. 6.2.6. 
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Choosing the point z to be z + , this implies 
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sm 

cos 



zi 
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Thus e.g. 
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-l^cos 
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L z () 



(E.9) 



provided the Bohr-Sommerfeld-Wilson quantisation condition holds, i.e. 



j\z(^- q hl-\h\U{z)^ =(2n + l)| 
where it is understood that (cf. Eq. (83)) 



n = 0,l,2,3,..., 



(E.10) 



-qh\ ~E-V{z±). 



Similarly under the same condition 
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We rewrite the quantisation condition in the present context and in view of 
the symmetry of the potential in the immediate vicinity of z + as half that of 
Eq. (jElO|) . i.e. 

1/2 



r z+ fx \ 1 7r 

J dz^-qhl--hiU(z)) ~q~. (E.ll 



The integral on the left can be approximated by 



^1 dzill 



■-0 



Z - Z 4 



1/2 



q . _! 2:+ - z 
2 Sm (2g//i 2 + )V2 



9 
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4 
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in agreement with the right hand side (the last step following from Eqs. (|119|) 
and (|120Jl ). We now see that Eqs. ()E.6|) and ()E.7|) assume a form as in our 
perturbation theory, i.e. they become 

tanj (g + 1)7 1 ~ Ttz, cot j (q + 1) T 1 ~ T^- (E.12) 



4 J ' 2j [ yi ' 4 J 27 

The equations corresponding to these in Ref. ^Hj are there described as 
WKB quantisation conditions of the double minimum potential. Since 

tan|(g+ 1)^| = for q = q = 3, 7, 11, . . . , 

cotj(g+ l)f} = for g = g = 1, 5, 9, . . . , (E.13) 

we can expand the left hand sides about these points and thus obtain 
27 CSOJ IT(hl)*ft(2z+)* i hl 4 



<?-<?o-T— = =f4\/- 77 77fj — e ~ 2 +z+ for g = 1,3,5,... 

7T7 V 2tt [i(g - 1)]! 

(E.14) 

in agreement with Eq. ()116j) . We note here incidentally that this agree- 
ment demonstrates the significance of the factor of 2 in the first equality of 
Eq. (jE.14|) which results from the factor of 2 in front of the sine in the WKB 
formula (EH). 



F Evaluation of WKB exponential with ellip- 
tic integral 

Here we evaluate the elliptic integral ^(0) of Eq. (|135a|) . 

F.l Parameters and their expansions 

The expression to be evaluated is 

(1 + u) 1/2 [E(k) — uK(k)\. (F.l) 



z 3G 2 

We have 



1—7/ , 

k 2 = , u = Gy/2q, g~2n+l. (F.2) 

1 + u 
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Hence we obtain for G close to zero: 
2 



and 



and 



k 



[l+u 



1 - 2Gy/2q + AG z q 



l + G^/2q~ 

[l + Gy/2i)W = l + G ) ll-\<?q 



, k*~l, 



1/2 



fc" 8 = 1 - k 2 = 2G^2q - 4G 2 q + • • • , k ~ 0. 
We now reexpress various quantities in terms of u. Thus 



k = 2u-2u 2 + ■■ 



and hence 



k' = v^(l - «) 1/2 = V^fl - I + • • • J . 



(F.3) 

(F.4) 
(F.5) 

(F.6) 
(F.7) 



The following expression appears frequently in the expansions of elliptic in- 
tegrals. Therefore it is convenient to deal with this here. We have 



k' V2«(l - f 



Hence 



ln(i|=ln 



1 + ^ 
2u\ 2 



Phl\ 2 



In 



:f.8) 



u 
2' 



;f.9) 



F.2 Evaluation of elliptic integrals 



Our next objective is the evaluation of the elliptic integrals E(k) and K(k) 
by expanding these in ascending powers of k' 2 , which is assumed to be small. 
We obtain the expansions from Ref. [2E] as 



E(k) = l + - 



k' 2 + - 

16 



mil)-- 

k' 12 



k' 4 + 



and 



k' 2 + .... 



(F.10) 



(F.11) 
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Consider first E(k): 



E{k) 



1 

1+ 2 



hi 



+ 



16 



In 



u 1 
h 2"2 
u 13 
2 ~ 12 



2«(1 - w) 
4w 2 (l -w) 2 + 



1 + ln 



u(l - u) + -u 2 (l - uf 



We can rewrite this as 



E(k) = 1 + In ( -^L ) M(1 - M) {4 + 3^(1 - u)} 



+— u{u - 1)(1 - «){24 - 3u(6u - 13)}. (F.12) 



Analogously we have 



uK(k) 
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+ --1 

2uJ 2 



2m(1 
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, «+ -«) } + 



l + (l-u)(- 



U 



in' 4 ^ 



-K(l-u) + 2u} + i u 2 [2 + (l-n)(n 



- 1 

2)] 



or 



uK(k) 



In 
In 



4 =^Hl-u) + 2} + in 2 [-n 2 + 3n] 



= ] - w) + 2} - -u\u - 3). 



(F.13) 
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With (|Fl2j) and (lF~T3j) we obtain 



E(k)-uK{k) 
4 

1 + ln 



- u){4 + 3«(1 - u)} - l -u{u{l - u) + 2} 
- 1)(1 - «)^{24 - 3n(6n - 13)} + -^{u - 3) 
1 + In (-^=) ^[(1 - M ){4 + 3u(l - u)} - 2{u(l - u) + 2}] 



+ -^—u(u - 1)(1 - u){24 - 3u(6u - 13)} + \u z (u - 3). (F.14) 
2.4.6 4 

Now consider the last line here without the factor 2.4.6 in the denominator 
and in the last step pick out the lowest order terms in u (i.e. up to and 
including u 2 )\ 

u(u - 1)(1 - u){24 - 3u(6u - 13)} + 12u 3 (u - 3) 
= -u(l - u) 2 {24 - 18u 2 + 39u} + 12u 3 (u - 3) 
= -3u(l - m) 2 {8 - 6m 2 + 13m} + 12u 3 (u - 3) 
~ -3w.8 + 6w 2 .8- 3.13m 2 

= -3u.8 + 9u 2 . (F.15) 
Now consider the bracket [...] in (|F.14j) . i.e. 

[(1 - u){4 + 3w(l - «)} - 2{u(l -u) + 2}] = u(3u 2 -4u- 3). (F.16) 
From (IF~T3l) with (lF~TIJ) and (lF~l5|) we now obtain 

E(k) - ~ 1 + In f i n 2 (3n 2 - 4u - 3) - -u + (F.17) 



F.3 Evaluation of 7 2 

We now return to Eq. (|F.1|) . i.e. 

2 



^(l + u) 1 / 8 ^*)-^)]. 
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Inserting here from the above expansions the contributions up to and includ- 
ing those of order w 2 , we obtain: 
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3G 2 



u 
2 
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In 



1 3 « 2 3n 



It 



2it/ 4 



(3 + Au - 3u 2 



(F.18) 



Remembering that u = G^/2q, this becomes 
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(F.19) 



Thus 
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(F.20) 
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